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PREFACE 


The  equations  governing  the  flow  of  a  conducting  gas 
interacting  with  a  magnetic  field  are  coupled  nonlinear 
partial  differential  equations  in  several  dependent  and 
independent  variables.  Thus,  exact  solutions  are  exceed¬ 
ingly  rare.  One  class  of  problems  which  can  be  solved 
exactly  involves  the  flow  of  an  infinitely  conducting  gas 
when  the  velocity  and  magnetic  fields  are  aligned  every¬ 
where.  A  particular  boundary  value  problem  involving  such 
a  flow  is  solved  in  this  Memorandum.  The  results  presented 
here  should  be  of  use  to  all  who  are.  working  in  the  area 
of  magnetohydrodynamics. 
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SUMMARY 


When  an  infinitely  conducting  fluid  flows  under  the 
influence  of  a  magnetic  field  in  such  a  way  that  the  velocity 
and  magnetic  fields  are  aligned,  then  the  governing  equations 
may  formally  be  reduced  to  those  of  ordinary  gas  dynamics. 

The  equations  are  still  nonlinear,  but  if  one  restricts 
attention  to  plane  flows, then  the  hodograph  technique  may 
be  employed.  The  equations  in  the  hodograph  tlane  are 
linear  ar  '  may  be  solved  exactly.  The  difficulty  then 
remaining  is  to  satisfy  the  boundary  conditions.  Only  for 
restricted  classes  of  problems  can  this  be  done.  This 
paper  treats  one  such  class  of  problems:  the  flow  of  a 
jet  of  gas  out  of  a  slit  in  a  rectangular  channel.  The 
method  used  is  a  modification  of  the  original  method  of 
Chaplygin  for  nonconducting  gas  jets. 
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THE  FLOW  OF  A  CONDUCTING  GAS  JET 
WITH  ALIGNED  MAGNETIC  FIELD 

1 .  INTRODUCTION 

In  this  Memorandum  we  treat  a  particular  boundary- 
value  problem  involving  the  flow  of  an  inf initely— conducting 
gas,  moving  in  such  a  manner  that  its  velocity  field  is 
everywhere  aligned  with  the  resultant  magnetic  field.  In 
such  a  situation  the  equations  governing  the  motion  can  be 
reduced  to  those  of  ordinary  gas  dynamics  [1-4,  6]. 

While  this  reduction  is  a  significant  one,  the  equations 
are  still  nonlinear  and  hence  intractable  for  most  boundary- 
value  problems.  If,  however,  we  limit  ourselves  to  plane 
flows,  it  is  possible  to  make  use  of  the  hodograph  technique. 
The  effect  of  the  hodograph  transformation  is  to  linearize 
the  basic  equations;  the  penalty  paid  for  this  simplification 
is  that  for  most  problems  the  boundary  conditions  in  the 
hodograph  plane  are  not  known.  These  are,  however,  classes 
of  problems  for  which  this  difficulty  can  be  overcome;  the 
outstanding  example  of  such  problems  is  the  efflux  of  gas 
from  a  slit,  the  jet  problem,  first  solved  by  Chaplygin. 

Chaplygin  solved  the  problem  of  jet  flow  from  a  slit 
in  a  half-plane,  for  which  the  velocity  infinitely  far 
upstream  must  be  zero.  Recently  this  solution  has  been 
extended  to  a  rectangular  channel  with  a  finite,  nonzero, 
upstream  velocity.  We  shall  be  concerned  with  extending 
this  solution  to  the  magnetohydrodynamics  case  when  the 
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gas  is  inf initely— conducting  and  the  velocity  and  magnetic 

fields  are  aligned.  The  velocity  at  any  point  v;ill  be 

2  2 

assumed  to  be  subcritical  (A^  +  M  <  1) . 

The  problem  can  be  reduced  to  the  solution  of  the 
equation  for  the  stream  function  in  the  velocity  plane 
under  appropriate  boundary  conditions.  The  solution  is 
obtained  by  separation  of  variables.  The  resulting  second- 
order  ordinary  differential  equation  for  the  velocity 
dependence  contains  the  hypergeometric  equation  as  the 
limiting  case  for  zero  magnetic  field.  One  of  the  two 
linearly  independent  solutions  of  the  equation  is  analytic 
for  velocities  close  to  zero;  the  series  expansion  of 
this  solution  about  zero  velocity  is  discussed  in  the 
Appendix. 


2. _ GOVERNING  EQUATIONS 

The  steady  flow  of  an  infinitely  conducting,  inviscid, 
perfect  gas  in  the  presence  of  a  magnetic  field  is  governed 
by  the  following  set  of  equations: 


(1) 

Continuity 

7-fpq) 

«  0, 

(2) 

Momentum 

pq-'q 

=  -  7p  +  i(7  X  B) 

(3) 

Energy 

q*7^) 

p  ■ 

•v 

o 

II 

(4) 

(5) 

Maxwell ' s 
Equations 

|  ~  ■  B 

x(q  x  B) 

Lx 

II  II 

o  o 

• 
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We  assume  B  is  parallel  to  q  and  write 

(6)  B  =  aPq  . 

(If  B  and  q  lie  in  the  same  plane, then  as  a  consequence  of 
Eq.  (5),  q  x  B  is  a  constant.  From  this  it  follows  that 
if  q  and  B  are  parallel  anywhere,  they  are  parallel  every¬ 
where.  Thus,  the  condition  expressed  by  Eq .  (6)  is  not  as 
restrictive  as  it  might  seem  at  first  glance.) 

Equations  (1)  and  (4)  lead  to 

(q  •  7)a  =  0 

i.e.,  that  a  is  a  constant  along  streamlines. 

Equation  (3)  is  satisfied  automatically  by  the  choice 
of  B.  The  momenti a  equation,  Eq.  (2),  may  be  written  as 

(7)  7($P  +  r  q  X  (7  /  q)  —  jjt  B  X  (7  x  B)  . 

•  4 

If  we  take  the  scalar  product  of  this  equation  with  q  we 
obtain 

(8>  |!+r^  =  Ho, 

where  Hn,  the  total  enthalpy,  is  constant  along  each 
streamline. 

From  this  point  on  we  shall  assume  that  each  of  the 
quantities  a,  HQ,  and  p/pY  =  f(S)  is  constant  throughout 
the  flow  field,  not  just  along  streamlines.  This  would  be 
the  case  if  all  the  streamlines  originated  in  some  region 

where  uniform  conditions  prevail.  In  the  particular  problem 


to  be  treated  in  this  paper,  such  a  situation  does  exist  and 
thus  the  indicated  assumption  is  not  actually  a  restriction. 

With  H0  constant  throughout  the  flow  field,  Eq.  (7) 
yields 

(9)  q  x  (7  x  q)  =  ~  B  x  (7  x  B)  . 

If  we  now  assume  that  the  flow  is  two-dimensional  (plane) 
and  write 


q  =  uT  +  v  j , 
then  Eq.  (9)  leads  to 

(10)  j|[v(l-Cp)]  -  |^[u(l-CP)], 

2 

where  C  =  a  /u  *  const,  (throughout  flow  field). 

For  plane  flow  the  continuity  equation,  Eq.  (1), 
bee  ,..es 

(11)  +  2|£vl  =  0. 

The  governing  equations  of  the  flow  have  been  reduced 
to  the  three  equations  (8),  (10),  and  (11).  Equation  (10), 
which  is  of  the  form  of  an  irrotationality  condition,  can  be 
identically  satisfied  by  introducing  a  potential  ce(x,y) 
defined  by 

lx  =  U(1“CP)  * 

(12)  1. 

,  vr;  =  v(l-Cp)  . 
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Similarly,  Eq.  (11)  can  be  satisfied  by  introducing  the 
usual  stream  function  Y(x,y)  such  that 


(13) 


far 

j 

1  BY 


> 


v  . 


The  three  governing  equations  may  now  be  combined  to  yield 
a  single  equation  in  either  Y  or  cp: 


(14)  Y2Y  +  2Y  Y  Y  +  Y2Y  -  (Y2+Y2) 
v  '  x  xx  x  y  xy  y  yy  v  x  y' 


( 1— CP^(M  — 1) 


(Y  +Y  )  =  0, 
v  xx  yy' 


2 

(15)  c  *t  +  2 co  cp  cp  +  cp  c?  ~  (rp  +co  ) 

v  '  xxx  xTy  xy  y  yy  v ^x  y' 


1+CP  (M  -1) 

L  M2 


(>:x+  yy>  “  °‘ 


(See  Imai  [1],  Kida  [2],  Seebass  [3,4].)  In  the  limit  of 
zero  magnetic  field,  C  -  0,  these  equations  reduce  to  the 
corresponding  equations  in  ordinary  gas  dynamics. 

A  principal  difficulty,  the  nonlinearity  of  Eqs.  (14) 
and  (15),  may  be  eliminated  by  interchanging  independent 
and  dependent  variables,  that  is,  by  transforming  from  the 
physical  plane  to  the  hodograph  plane.  Proceeding  in  the 
usual  manner,  we  obtain  the  following  equations  for  Y  in 
terms  of  the  new  independent  variables  q  and  9,  polar 
coordinates  in  the  hodograph  plane: 


(16)  q^(l~CP)Yqq+ 


(1-CP)2(1+M2)  +  CP  M^[3— CP+y(CP-1)  ] 
1  -  CP  ( 1  -  M2) 


f  (1-M2)  [l-PC(l-M2)]Ye9  =  0. 
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An  equivalent  equation  can  be  written  down  for  x,  but 
since  this  will  not  actually  be  required, we  shall  not  set 
it  down.  Once  Y(q,9)  is  determined  from  Eq .  (16),  (q,9) 

may  be  determined  from  the  hodograph  relations,  which  in 
this  <~ase  take  the  form 


(17) 


q(’-cp)2,.  T  ( 

p  l-CP(l-M^)  q 


< 


co_  = 


-  4  (l-cnn-iC) 


F  2 

2'  ' (1-CP)  +  CPMZ 
l-CP(l-M^) 


Y, 


Is 


2 

(M  and  p  can  be  expressed  in  terms  of  q  by  use  of 
Bernoulli's  equation,  Eq.  (8).) 

Once  Y(q,9)  and  (q,9)  are  known,  the  coordinates  in 
the  physical  plane  can  be  obtained  from  the  expressions 


1  Pr 

dx  =  q'ci-cp)  (cos  9  dcv  ~  ^r(1_cp)  sin  9  df)  ' 


(18)  1 


1  p 

dy  =  q(l-^p)'^  -^(l-CP)cos  9  dY  +  sin  9  d  )  . 

vy 

3_. _ FORMULATION  ,^ND  SOLUTION 


The  particular  problem  to  be  solved  in  this  Memorandum 
is  the  outflow  of  an  infinitely  conducting  gas  from  p  slit 
in  a  rectangular  vessel.  This  problem  has  been  solved  for 
the  nonmagnetic  case  by  Fal'kovich  [5].  The  configuration 
is  shown  in  Fig.  1.  The  channel  is  assumed  to  be  of  width 
H,  the.  opening  of  width  h.  Infinitely  far  upstream  the 


and  is  aligned  along  a 


flow  is  uniform  with  velocity  Vq, 
uniform  magnetic  field  Bq .  Since  the  flow  velocity  and 
magnetic  field  are  aligned  at  upstream  infinity,  they  remain 
aligned  throughout  Lae  flow  field,  as  indicated  earlier. 
Also,  a,  H  ,  and  S  may  be  taken  constant  throughout,  since 
all  the  streamlines  originate  in  a  region  of  uniform  con¬ 
ditions.  Thus,  the  conditions  that  are  needed  for  the 
theory  presented  in  Sec.  2  to  apply  are  satisfied. 

Chaplygin's  method  for  two-dimensional  gas  jets  fails 
when  the  flow  in  the  jet  is  partly  subsonic  and  part'  r 

supersonic  (see  [7]).  Since  the  hypercritical  transition 
2  2 

(A  4-  M  -  1)  has  the  same  essential  character  as  the  sonic 

transition, one  would  expect  that  in  the  present  problem  the 

restriction  to  be  imposed  is  that  the  flow  remains  sub— 

2  2 

critical  everywhere,  that  is,  A  +  M  <  1.  We  shall  assume 
this  to  be  true  ir.  what  follows. 

The  jet  is  bounded  by  two  free  boundaries,  C 'D  and 
CD,  on  which  the  pressure  is  constant.  The  stream  function 
Y  is  assigned  the  value  along  the  upper  boundary  of  the 
flow,  A3'C'D,  and  the  value  —  along  the  lower  boundary, 
ABCD,  the  total  mass  flow  rate  then  being  Q.  The  final 
velocity  far  downstream  on  the  free  boundaries  is  v^. 

The  physical  flow  plane  represented  in  Fig.  1  can  be 
mapped  onto  tne  hodograph  plane.  In  the  hodograph  plane 
we  introduce  polar  coordinates  9  and  t,  where  t  =  q  /q  max* 
The  limiring  values  of  t  far  ”pstream  and  downstream  are 
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then  given  by  tq  =  v2/ q2max  and  ^  =  vj/q2max.  The 
hodograph  representation  of  the  flow  is  shown  in  Fig.  2, 
in  which  corresponding  points  are  designated  by  the  same 
letters.  CDC 1  is  a  semicircle  of  radius  EAE1  a  semi¬ 
circle  of  radius  t  .  The  values  that  Y  must  take  on  the 

o 

various  boundaries  in  the  hodograph  plane  are  given  below: 


(19.1)  Y  =  —  -jQ  when  T  =  T^,0<9^-jTr  DC 

(19.2)  Y  =  -  -JjQ  when  £  t  >  0,  0  =  ^rr  CB 

(19.3)  Y  =  -  when  0£t<tq,  9=40  BA 

(19.4)  y  =  +  When  t  =  Tp  -  ^tt  £  9  <  0  DC* 

( 19 5)  Y  =  +  when  T^^.T>0^5=--2Tr  C  *  B  1 

(19.6)  Y  =  +  when  O^t^Tq,  9=-0  B'A 


The  problem  now  is  to  find  a  solution  of  the  basic 
equation  for  Y,  Eq.  (16),  satisfying  the  boundary  conditions 
(19.1)  -  (19.6)  in  the  hodograph  plane.  Before  attempting 
this  we  must  first  transform  Eq.  (16)  from  independent 
variables  q,  9  to  the  new  variables  t,9«  We  also  note  chat 
the  Alfven  number,  defined  as 


(20)  A 


flow  speed 
Alfve;n  speed 


is  related  to  the  constant  C  as  follows 
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(21) 


CP  = 


1  _  azp 

7  =  ~ 


Then,  introducing  the  stagnation  Alfven  number  Aq , 


(22) 


*  _  1  fh 

0  ‘  °  K  ’ 


we  can  write  (assuming  a  perfect  gas) 


(23) 


A2  % 


-  (1-T) 


(24) 


-j-a-T) 

Ao 


Also, 


(25) 


p  =  po(1~t) 


M2  2  T 

M  -  t=f  • 


Now  making  the  change  of  variables  from  q,  9  to  t,  h  and 
introducing  relations  (24)  and  (25),  Eq.  (15)  transforms 


into  the  equation 


r  1 1  r  2-y  ~ 

4^2(1-t)  L1  ~  A^2(1-t)  L1  +  A~2(1-t) ^(^-l)  *f TT 
fr  nr  i 

(26)  +  4t  |1  -  Aq2(1-t)Y=T  i  1-T )  ( 1-A~2 ( 1-T ) ^=T)  + 

^  2  r~  n 

-  2a0-2t2(1-t)^  1=2  ^  A0-2(1_T)Fll  I  Tt 

r  (v  u  2-y  i,  H 


+  UL-^=Tt  1  +  A02(1-t)y-I(^3  t-1)  y,9  =  0 
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We  must  now  seek  a  solution  of  Eq.  (26)  satisfying  the 
boundary  conditions  (19 . 1)— (19 .6) .  We  shall  do  this  by 
finding  separate  solutions  in  these  subareas,  lauelled  (1), 
(2),  (3)  in  Fig.  2. 

For  regions  (1)  and  (2)  we  look  for  separable  solutions 
of  the  form 


(27) 

(t, 9)  =  -  ^ 

4 

cr 

2  anGn(T)sin 
n=l 

2n9, 

(28) 

¥^2^(t,9)  =  4  ^ 

■  + 

05 

2  a  G  (r)sin 
n=l  n  n 

2n9 . 

Substituting  either  one  of  these  expressions  into  Eq. 
(26),  we  obtain  the  following  equation  for  Gn(t) : 


T  ( 1“ T )  l-V(l-r) 


1  +  A" 


+  1-A  2(1-t)v 


r 


(29) 


J  l_ 

2— Y 


l^o  d-T)Y  + 


-2Ao2  U  ,  +  ^  Ao2(1-t)^  l 


JCv-1) 


"2  0  t-1  i  1  +  A02(1-t)^~^ 


\Y_1 


Gn  =  0. 


We  shall  consider  this  equation  in  greater  detail  later. 
For  the  present  let  the  two  linearly  independent  solutions 
of  th  equation  be  denoted  by  Gn(T)  and  H^t)  .  It  will  be 
shown  later  that  only  one  of  the  solutions  remains  bounded 
at  t  =0;  let  Gn(r)  be  this  solution. 
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The  expression  assumed  for  Y^^(t,9)  satisfies  boundary 
condition  (19.3)  and  boundary  condition  (19.2)  on  sector 

(9) 

BE,  while  the  expression  assumed  for  t'v*'/(t,9)  satisfies 
boundary  condition  (19.6)  and  boundary  condition  (19.5)  on 
sector  B 'E ' . 

In  region  (3),  which  represents  the  annular  region 
CDC'E'AEC,  we  look  for  a  solution  or  the  form 

(30)  t(J)(T,9)  -  -  §9  +  2  [AnGn(T)  +  BnHn(r)]sin  2n9 . 

n=l 

C3) 

This  expression  for  Yv  '  satisfies  Eq.  (26)  identically; 
it  also  satisfies  boundary  conditions  (19.2)  and  (19.5)  in 
sections  CE  and  C'E1,  respectively.  To  satisfy  boundary 
condition  (19.1)  and  boundary  condition  (19.4)  on  arc  CDC 1 
(on  which  t  «  t^)  we  require 

°° 

_  +  s  [AnGn(Tj)  +  BnHn(T1)]  sin  2n»,  0  <  9  £  in  , 

n=  1  ^ 

=  _  g9  +  z  tAnGn<Tl>  +  BnHn(Tl)]  sin  2n9,  -  ^  9  <  0  , 

1  n  n=l 

or  equivalently, 

(31)  2e  +  =  2  !\Gn(T1)  +  BnHn(T1))  sin  2n8, 

n=i 

where  the  upper  sign  is  taken  when  9  >  0,  and  the  lower  one 
when  9  <  0 . 
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Since  the  left-hand  side  of  Eq.  (31)  is  an  antisymmetric 
function  of  9  in  the  interval  —  ^  ®  ^  -jtr,  it  can  be 

expanded  as  a  Fourier  sine  series 


(32)  S8  +  tQ  -  -  3  “  H  sin  2n9>  '  7"  ^  !  <  I’ 

n=l 


Equating  the  right-hand  sides  of  Eqs.  (31)  and  (32)  we 
obtain 


(33) 


A  G  (T.)  +  B  H  (t.)  = 


_  2_ 

nn 


This  completes  the  satisfaction  of  boundary  conditions  (19.1)- 
(19.6)  . 

There  are  additional  boundary  conditions  which  must  be 
satisfied.  These  result  from  the  requirement  that  the 
expression  for  ¥^(t,9)  be  the  analytic  continuation  of 

and  Y^)  into  region  (3).  This  leads  to  the  conditions 


(34) 


'?(3)C„,9)  - 


f(3)(T0,e) 


f(1)(ve)’  H(\<9>  -  I?(1)(to’9>'  0  i  9  £ 


('2')  t^v(3)  p^v(2)  ^ 

tw<t0,9),  It  (T„,»)-It  <v')*  9  1° 


Substituting  the  expressions  for  ^ y(^) ^ ^(3)  ^nto  £qS> 
(34)  we  obtain 


r ' 


(35) 


-  §9  +  2  [A  G  (t  )  +  B  H  (t  )J  sin  2n9 

tt  ,  1  n  r.  o  n  nv  o'  1 

n=  l 


co 

7  +  x  a  G  (t  )sin  2nQ,  0  <  9  £ 
n=l  n  n 


2  [A  G'(t  )  +  B  H'(t  )]  sin  2n9 
n_-^  1  n  nx  o7  n  n v  o'  J 


=  2  anGn<To)  sin  2nP' 

n=l  n  n  ° 


0  1  1  f*> 


(36) 


GO 

-  5s  +  2  [A  G  (t  )  +  B  H  (t  )  ]  sin  2n<> 


< 


tt  ,  ‘  n  n'  o'  n  n '  o' 

n=i 


+  n^anGn^To^  sin  2n°* 


2  [A  G'(t  )  +  B  H'(t  )]  sin  2n« 
-i  1  n  n  ^  o'  n  n v  o'  J 
n=l 


f  ^  1  0, 


=  2  a  G  v  )  sm  2nQ , 

„  n  n  nN  o'  * 

n=  I 


-  fr  <  9  £  0 


w 


Combining  the  first  equations  from  Eqs.  (35)  and  (36),  we 


can  write 
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where  the  upper  sign  is  taken  when  9  >  0,  the  lower  one 
when  <  0.  As  before,  the  left-hand  side  can  be  written 
as  a  Fourier  sine  series  (Eq.  (32));,  Eq.  (37)  then  reduces 
to 


(38) 


<An- 


a  ) 
n7 


G  (t  ) 
nv  o7 


B  H  (t  )  = 
n  nv  o' 


_ 

~n 


The  second  equations  of  Eqs.  (35)  and  (36)  yield  the 
single  relation 


(39)  (VVW  +  BnH^(To)  =  0. 


Equations  (33),  (38),  and  (39)  determine  the  three 

sets  of  unknown  coefficients:  a  ,  A  ,  and  B  .  Solving 

n  n*  n  ® 

these  equations, we  find 


a  = 
n 


(40)  < 


A  = 
n 


B  = 


-  JL  _  1 

7111  Gn(TP 


nn  Gn(Tl^ 


2_  Gn(To> 

"n  W  ' 


G'(t  )H  (t,)  -  G  (t.)H'(t  ) 
^  +  _n__o_  nv  l7  nv  l7  nv  o7 


r 

i  + 

L 


G 1 (r  ) H  (t,) 
n  o7  nv  l7 


W 


where  VL(t  )  is  the  Wronskion  of  the  solutions  G (t),H  (t) 
no  n  n 

evaluated  at  t  ,that  is 
o  ’ 
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(41)  W  (t  )  «  G  (t  )H'(t  )  -  G'(t  )H  (t  ) 
v  '  nv  o'  nv  o'  nv  o'  nv  o'  n  o' 


If  we  now  substitute  these  results  into  the  expressions 
for  Y^^,  Y^3\  the  solution  takes  the  form 


(!)  n  0  “  1 

(42)  T  (T,8)  .-§-2  a 

n=l  n(T,) 


Gn(To>Hn(Tl)  "  Gn(Tl)Hn(To> 

1  +  n  °  n  /,  ,  G  (t)  sin  2n9, 

wnk  o'  J 


(43)  Y(2)(t,9)  = 


uu 

-  £  2 


n  n=l  nVTl> 


!  ,  GA<To>Hn<V  "  Gn^l)H^To> 

i  + - 7~, — 7~Z — V - 


n  O' 


G  (t)  sin  2n9, 

4  1 


(44)  Y(3)(t,9) 


_2e-9  £  L  1  _  Jl  +  Gn(T°>Hn[li) 

r  "  L  w n(T„) 


w^;  G-(T) 


j  ^4^7  I7  VT>  |  sin  . 
L  n  0  J 


Equation  (44)  may  also  be  written 


(45) 


{j  ’,(3)(t,9)  = 


-  9  -  Z 


X„(t) 


sin  2n9, 


where 


17- 


(46) 


G  (T) 


n(T)  -  WnrpT 


1  VT>Gn<Tl>  -  Hn<Tl>Gn<T> 


1'  "n 


tiT 


Now  that  the  solution  in  the  hodograph  plana  has  been 
obtained, it  is  necessary  to  transform  back  to  the  physical 
plane.  If  we  introduce  the  hodograph  equations,  Eqs.  (17), 
into  the  transformation  relations,  Eqs.  (18),  we  obtain 


dx  mk 


~  P  C  B  p 

-  ^  sin  8  '?0  +  - ^  p— q(l— CP) 

-  1— A 


1  -  CPU  +  §  g) 


-1 


de 


pTsin  8  Yq  +  [q(l-CP)  §q(  p~  )  V 


dq 


dy  =  -l 


1— A 


"in_?  q(l-CP)  (p|)Y9  +  cos  9  Yf 


dq 


jr  q(l-CP)2H  -  CP ( 1  +|  fj|)]  1  yq  +  ^os  9Y( 


1— A 


de 


> 


Integrating  these  expressions  results  in 
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y  =  q 


Z77  t— CP ( 1+  ^  -r~)  I  !  v  sin  9d9  +  ■—  cos  fid 

1— A  i")  w  P  clcl  -*  L  a  r  h 


+  y0(q)  • 

Making  the  change  to  the  variable  t  and  introducing  the 
relations  Eqs.  (21)— (25),  we  finally  obtain  the  following 
for  the  coordinates  in  the  physical  plane: 


9 

_fyc  sin  9d9  +  2t[1  -  A~2(1-t)^] 


fl-23T(l-T)  1)]  1  rT  sin  QdQ 

T 

a 

+  /  l'a  r-os  «dQ  +  yo(r)  , 

where  for  convenience  we  have  introduced  0  =  l/(y—  D  . 

The  y-coordinate  in  region  (3)  can  be  found  by  sub¬ 
stituting  the  expression  for  Eq.  (45),  into  Eq.  (48). 
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We  obtain 


09)  y»)  =  -  J  -  S-i?ST±iliJ 


-]>. 


where 


(50)  R(t)  =  [L-A02(l-T)p)[l-Ao2(l-T)Pfl-2f?T(l-r)-1l)''1 


j.  a  -L  v  v  r  sin(2n+l)  9  ,  sin(2n- 
n^Xn(T)  1 "  2’n+l  +  — 7-Sri- 


(Since  y  =  0  when  9=0,  y  (t)  has  been  set  equal  to  zero 
in  Eq .  (49)  .) 

Along  the  outer  boundaries  of  the  jet  CD  and  C'D  the 
velocity  is  consta  %  equal  to  Tj_.  Setting  r  =  rL  in  Eq. 
(49),  we  obtain 

('3')  Q(l— 1 Ti)~^  *  X  (t,) 

(51)  y(3)(Tl,9)  =  -  __± - Jt  R(t  )  2  _n_J_ 

M1  1  n=l 


[Sin(2n-1)°  _  sm(2n+l)9, 
1  in-1  Zn+T 1 


+  sin9  +  2  Xn(TL)  [ 


sin(2rH-l)9  ,  sin(2n-l)9 
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From  che  expression  for  X  (t),  Eq.  (46) ,  we  see  that 
Xn(T^)  =  1.  The  Last  sura  in  the  above  equation  can  then 
be  written 

v  r  sin(2n+l)  9  ,  sin(2n-l)9,  n  .  0  v  sin(2n+l)9 

nfp  2tvfl  L  +  Zn-1  1  =  sln  9  T  2nlx  ~inn  ~ 

The  right-hand  side  is  the  Fourier  sine  series  expansion 
of  the  function  —  sin  9  +  ^  (with  +  for  9  >  0,  —  for  9  <  0) . 
Hence  we  can  write 


(52) 


2 

ri=l 


rsinI2n 20-1)9, 
1  2n+T  ^  2thT j 


a  ,  n 
sin  9  +  tj ) 


+  for  6  >  0, 
-  for  9  <  0 


Substituting  this  expression  in  Eq.  (51)  yields 


(3)  Q(l~Ti) 

(53)  y(3)(Tr9)  ^ 


-o 


VlR(Tl)  s 


X  (T,) 

nv  1' 


n-1  " 


r sin(2n-l) 0  _  sin(2n+l) 9 , 
*■  2n— 1  2n+l  * 


±  -!  '<  ,  0  £  0  • 


When  9  -  —  •»,  is  equal  to  half  the  width  of  the  opening 

of  the  jet,  y,  and  thus 


(3)  r  h  Q(l— I 

y  )(’-1,  -  j)  -  7  f  iR(Ti>  n_1  n 


?  Xn(Tl>  sin  T 


L 


T5=r 


r\ 


sin  Qa filn 

- 2TSI  +  7>* 


or, 


(54)  h  =  MizIll 


-0 


-  (-i)nx;^Ti> 

^TlR<Tl)  s,  --T 


TT 


n=l  4n  -1 

l 

By  continuity  of  mass,  we  have 


P1  1 

(55)  Q  =  mass  flow  =  q-^h  , 

ro 


where  h  is  the  width  of  the  jet  at  x  =  +  ®.  Using  Eq.  (25), 
we  may  rewrite  (55)  as 


Q  =  h  * 

Using  this  expression  to  eliminate  Q  in  Eq.  (54),  we  obtain 


1 


E  " 


h 

r 

h 


i-jf 


t1R(t1) 


n'l  — ^cr~ 


where  k  is  the  coefficient  of  contraction  of  the  jet. 
Differentiating  Eq.  (46),  we  find 
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(56)  Xn(T1)  = 


G  (ti) 
n v  1' 


1  c7?p- 


VT1>  Gn<To> 

^^7 


Introducing  this  into  the  equation  for  1/k  above,  we  obtain 


(57) 


=  1 


8 

TT 


n 


1R<T1' 


z  L=n 

n=l  tXl 


°n<Tl>  VTl>  Gn<To^ 

Gn^TP  WnCTo*'  Gr/TP  j 


For  an  infinitely  broad  vessel  (H  -  !X)^  tq  must  be 
zero  and  the  expression  above  becomes  (denoting  the  value 
of  k  in  this  case  by  kOT) . 


1 


1  ~  ~  TiR(t-,) 


TT 


CC 

2 


n=l 


c-l)n 

"’T~~ 
4n  -1 


ten(Ti)  wn(TX)  Gn(0)  I 


This  should  be  compared  with  the  Chaplygin  formula  for  the 
nonmagnetic  case 

f  1  \  =  ,  _  8  T  2  (~P n 

Kj  nonmagnetic  -  1  ^  4n^l1  Z^T 

where 


Zn(r)  =  Tn  F(an,bn,  2n+l ;  t)  . 


APPENDIX 


The  equation  obtained  after  separating  variables  was 
shown  to  be  Eq .  (29),  that  is, 


I  2— v 

t2(1-t)  [l-A^d-r)^]  [1  +  a;2(1-t)^=T(^  T-l)]Gn(T) 


T  ! 


Y 


+  T  <  [1  -  A02(1-T)^1C]  [l-T-A^Cl-T)^1  +  Wp’ 


o 


(A— 1) 


L 


-  2A^2t2(I-t)^T[>^ 


1 


o 


1  ■n  —  - y 

(v-d7  7=1  ° 


1  A  2d-T)^] 


2-y 

+  n2  [XL*  t-1]  [1  +  A^CI-t)^^  t  1)  ] 2  Gn(r)  =  0. 


The  solution  we  require  is  the  one  which  remains  bounded 
when  t  =  0.  In  this  Appendix  we  shall  present  a  brief 
discussion  of  Eq.  (A— 1) .  Complete  details  will  be  given  in 
a  subsequent  paper.  First,  we  note  that  when  the  magnetic 

_9 

field  vanishes,  i.e.,  A  =0,  the  equation  reduces  to 
(A— 2)  t2(1-t)g"  +  t [1  +  (0-1) t JG^  +  n2 [ (1  +  28>T-l]Gn  =  0. 


This  is  the  usual  equation  which  arises  in  the  study  of  the 
hodograph  technique  [7].  The  solution  of  this  equation 
which  remains  bounded  at  t  =  0  is  given  by 


-24- 


(A-3)  Gn(T)  =  ^  F(an>bn'  2n+1’  T>  > 

where  a  ,  b  are  determined  from 
n  n 


(A-4)  an  +  bn  =  2n  -  0  , 


anbn  =  -Pn(2n+1)  > 

and  F(a  ,b  ,  2n+l;r)  is  the  hypergeometric  function,  whose 

1 1  Li 

power  series  . 


(A— 5)  F(an,bn,2n+1;T) 


2 

m=0 


r(an'fm)r(bn+m) 


is  uniformly  convergent  in  the  domain  jrj  <  1. 

When  n  =1,  Eq.  (A-l)  may  be  solved  exactly  (see 
[4]).  For  n  ^  1,  the  primary  difficulty  in  solving  the 
equation  resides  in  the  fact  that,  for  general  v,  the 
coefficients  involve  nonintegral  powers  of  t.  Note,  how¬ 
ever,  that  if 

(A— 6)  — ^  =  k,  or  v  =  1  +  p 


where  k  is  a  positive  integer,  then  the  coefficients  are 
polynomials  in  t.  In  this  case  the  equation  becomes 
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T2(L-')  [l-A0?(l-T)k][l+A02n-')k  L((2k+l)r-l)  ]Gn 
+  t  f [ 1-A^2 ( 1-") k ] t 1-t-a;2 ( L-t) k+1+kT ) 

(A— 7) 

-  2A;2T2(l-T)k-’-[(k-2k2)  -  kA~2(L-'r)k] 

+  n2 [ (2k+l)T— 1  ]  [  1+A~2(1— T>k— 1  ( (2k+l)  1)  ] 2  GR  =  0. 

The  finite  singularities  of  this  equation  are  given  by 


t  =  0, 


(l-')k_1[l  -  (2k+l) r J  =  A2. 

Since  k  is  a  positive  integer,  the  last  expression  has 

exactly  k  roots.  Thus  F.q.  (A— 7)  has  (3+k)  finite  singularities. 

In  the  special  case  k  =  1  (v  =  2)  the  four  finite  singularities 

are 

t  =  0, 

T  =  1, 

(A— 9)  2 

T  =  1  —  Aq  > 

'  =  3<^o>  • 

One  can  show  that  t  =  ®  is  also  a  singular  point. 
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Although  t  «  0  is  a  singular  point  of  the  general 

equation,  Eq.  (29)  or  (4—1),  this  equation,  like  the  hyper— 

—2 

geometric  equation  to  which  it  reduces  when  Aq  =0,  has 

one  regular  solution  near  t  =  0.  It  is  this  solution  which 

has  been  called  G  (t) .  To  find  this  regular  solution  we 

write  G  (t)  as 
n v  ' 


Gn<  >  =  2  ankT  ’ 

and  substitute  in  Eq.  (A— 1) .  At  the  same  time  the  co¬ 
efficients  of  the  equation  are  also  expanded  in  series  in 
t.  The  following  recurrence  relation  for  the  coefficients 

a  ,  is  then  obtained  : 
r.K. 


(l-A;2)2(k2-n2)ank  _  [  (k-1)  (k-2)  +  Aj  0-2  +  3Aq2  -  (y+1)A^*1  (k-1) 


(A— 10)  -n2r£!  +  -ij.  Aq2[-6v  +  Ao0V,-l))]J  an  k_1  +  n 


.-2 


2  a;2 

<  -i)2 


[8v2  -  7y  -  Aq2(12y2  -  U(  +  l)]an  k_2  , 


+  A 


-2 

'o 


s=0 
1= k 


s=k— 1 
't+s=k 


C,  s(s-l)a  - 
l  v  /  ns 


s=0 
_=k 

> 

7^ 
s=k— 2 
i+s=k 


sa 


ns 


+ 


n 


s=0 

h—k 

s=k-3 
1+ s=k 


cf  a  e 

ns 


■  o,  k  £  0, 


where 
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(A-ll-a) 


(A-ll-b)  C? 

u- 


(A— 11— c)  C  l 
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and  C\  ^  '  is  defined  in  the  usual  way 


=£)■ 


n, 


(  n— m) 


,’m! 


According  to  (A-10)  each  a  ^  is  given  in  terms  of  all 
previous  a  whereas  in  the  nonmagnetic  case  the  recurrence 
relation  involves  a^  and  only.  Thu^,  it  seems 

impossible  to  write  down  the  various  a  ^  explicitly.  How¬ 
ever,  we  do  note  that  the  first  nonzero  coefficient  occurs 
when  k  =  n.  It  is  then  more  convenient  to  write  our 
solution  as  follows 


(A— 12)  G  (t)  =  Z  a  Tk  =  Z  a  .  Tk  -  Tn  Z  a  Tk-n 

n  k=0  nk  k=n  nk  k=n  nk 


n  v 
=  t  2  a 


•  P  -  ^ 


p=0  n,p+n 


2  a 
5=0  nP 


The  recurrence  relation  for  a  ,  obtained  directly  from 
(A— 10) ,  is  found  to  be 


(l-A02)2p(p+2n)anp  -  U(p+n-l)  (p-l  n-2)+  ^(v-2  +  3Ao2-(y+1)A^ 


-  "2(£l  +  yTT  a/^yV^-^JVp-I  +  n 


,-2 

2  Ao 
(  Y-l) ' 


[8Y2-?Y-Ao2(12Y2-UY+l)]an(p_2 


)  (p+n— 1) 
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(A — 13) 


—  s=0 

+  Ao2 

t=p+n 

2  cj 

s(s— 1) a  — 

v  y  n , s— n 

s=0 

t~l 

s=p+n— 1 
s+-l=p+n 

s=0  ~ 

t»p+n 

-t=p+n 

2 

2 

C,  s  a^ 
v  n,  s 

+  n2  Z  C?  3  ,  „ 

— n  —  -t  n ,  s— n 

t»2 

t=3 

s=p+n- 

2 

s=p+n— 3 

s+t=p+n 

s+'t=p+n  - 

The  first  fev  terms  in  the  series  for  Gn(*r)  may  now  be 
calculated  from  (A— 13) .  We  find 


G  (t)  =  a  t°4  1  + 

/  nr\  \ 


'n'  '  no  '  '  (l-A~2)Z(l+2n) 


[  n  (  n-1)  +Bn-n2  D+n  (  n- 1 ) 


L 


(1— A~2) 222(l-fn) 


{[n(n+l)+B(n+l)-n2D+A^fn(n+l)  ] 


(A— 14) 


[n(n— 1)+Bn— n2D+n(n-l)  A^  J 

0-A~2)2a+2n) 


i—2 


A72n 


0  n(n-l)A^  ..  ..  ..  .. 

-  [n2E  +  - {y+A  2(1-3y)}  -  -2—  2(y— 3+A  ^(l 

(Y-l)2  0  (y-1)  2  0 


-3y)) 
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where 


B  =  T^T  [V-2  +  3  A  2  -  (v-H)A^i, 


(A-15)  B  =  Aq2(-6v  +  A02(5y-1))  , 


_0 
A  - 

E  =  — 2 — 
(Y-i) 


[Sy""  -  7y  -  Ao^(12y^  -  lly  +  1)  ] 


_2 

If  A  “  is  set  eoual  t  zero  in  (A-14)  the  series  reduces 


(A-16)  Gn(T)  ==  ano  rn  F(a  ,b  ,2i»fl-T), 


where  a„  and  are  the  same  quantities  defined  in  (A-4) . 

In  addition  to  Gn(r) ,  it  is  necessary  that  we  know 
W  (t) ,  the  Wronskian  of  the  two  linearly  independent 
solutions  of  (A-l) .  This  is  given  by 


(A-I7)  Wn(r)  =  WQ  iexp  -  jtq  44  dr  , 


where 


2-y 

Y^I  y+1 


a(T)  -  T(I-T)[l-A^(l-T)y-i]  [1  +  Aq  ( 1— t )  (S  t-1)]. 


b(T)  -  [l-A^l-T^j  [l-T-A2(l-T)V-'i  f^J 


-  2A;2  t2(1-t) 


2-y 

V-'I  r  V-3 


2  v-T  o 


4  a;2(i-t)^ 


( y — I)  "  T  "  w 

Equations  (A--14)  and  (A— 17)  are  sufficient  to  enable 


one  to  calculate  k,  the  coefficient  of  contraction  of  the  jet, 
using  Eq.  (37). 
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